Superconductivity in a single Ceo transistor 
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Single molecule transistors (SMTs) are currently attracting enormous atten- 
tion as possible quantum information processing devices [H E], SI SI [5]. An 
intrinsic limitation to the prospects of these however is associated to the pres- 
ence of a small number of quantized conductance channels, each channel having 
a high access resistance of at best i?x/2 = /i/2e^=12.9 kf^. When the contacting 
leads become superconducting, these correlations can extend throughout the 
whole system by the proximity effect and all contact resistance disappears. This 
not only lifts the resistive limitation of normal state contacts, but further paves 
a new way to probe electron transport through a single molecule. In this work, 
we demonstrate the realization of superconducting SMTs involving a single Ceo 
fullerene molecule. The last few years have seen gate-controlled Josephson su- 
percurrents induced in the family of low dimensional carbon structures such as 
flakes of two-dimensional graphene [6] and portions of one-dimensional carbon 
nanotubes [7j. The present study involving a full zero-dimensionnal fullerene 
completes the picture. 

The rapid experimental progress in SMTs has revolutionized our understanding of single- 
spin magnetism [31 lU [5] , nano-electromechanics [8[ |9] and Kondo correlations [TOl [HI [121 [13] • 
However, no experimental realization of superconductivity in a SMT has been reported so 
far. Superconductivity is a central paradigm to electronic phase coherence phenomena and 
its effects on SMTs transport properties have been extensively investigated by theory [HI [151 
[T6] . Josephson weak links relying on bottom-up quantum dot (QD) structures could recently 
be experimentally realized using semiconducting nanowires [17] as well as carbon nanotubes 
[71 [IHl [H]. Superconductivity in a SMT was yet to be observed as it offers new insights 
into the couplings involved in single molecule junctions. Superconducting correlations have 
been observed in ungated junctions involving the metallofuUerene Gd@C82 [20j but the 
absence of supercurrents have been attributed to the destructive effect of the local magnetic 
moment carried by the Gd atom on the superconducting phase coherence. In this Letter, we 
demonstrate the first experimental realization of superconducting single fullerene molecular 
transistors, in a broad range of coupling strengths to the contacts, covering both the weakly 
coupled sequential tunneling regime, as well as the resonant tunneling regime. 

SMTs are prepared following the electromigration technique [21j, further optimized as 



2 



presented in refs. (iQl [HI [121 [13]. By extending these advances to measurements at dilution 
refrigeration temperatures, precise electronic spectroscopies and measurements with super- 
conducting electrodes can be undertaken. The superconducting contacts to the molecular 
QD are either made of aluminum, or of gold in which superconductivity is induced by the 
close vicinity of an aluminum capping layer (see Methods for details). After electromigration 
performed at 4 K, stable devices exhibiting QD behavior are further cooled and measured at 
40 mK, with an effective electron temperature of 80 mK. We show here the data on 4 out of 
11 samples that were measured at low temperature and which exhibited Coulomb blockade 
and/or gate-modulated Kondo-like features. 

The superconducting weak coupling regime is illustrated in Figure [T}3, in a device (sample 
A) displaying two Coulomb diamonds (CDs). In the absence of other conducting regions 
near the Fermi level over the experimentally accessible gate voltage window (V^ = -4.5 to 
+4 V), an addition energy Eadd > 250 meV can be inferred. This brings clear evidence for 
single molecular quantum dot behavior because a metallic nanoparticle would lead to a much 
lower addition energy, typically not exceeding some tens of meV [22j. At base temperature 
and zero magnetic field, the conducting regions of device A do not intersect the Fermi level: 
they are separated by a spectroscopic source-drain voltage gap Q of about 680 /iV near the 
degeneracy point. 

This gap, refiecting the quasiparticule spectrum of the contacting electrodes, is a typical 
feature of a nanostructure weakly coupled to superconducting electrodes [23j. In such de- 
vices, the current is dominated by incoherent sequential tunneling of quasi-particles of energy 
|E| > Aq (where Aq = 177 /xeV is the bulk superconducting gap of aluminum). The extent 
Q of the non-conducting region is very close to 4Ao/e, as expected in an S-I-S tunneling 
structure (Figure [l]3,d). The source and drain contacts in a SMT usually have a strong geo- 
metrical asymmetry in electromigration devices (Figure [l^). In nanostructured aluminum, 
superconductivity can stabilize up to much higher fields than in bulk material, exceeding 
3 T in nanoscaled aluminum particles [22j; one therefore may expect very diflFerent critical 
fields in each contact. This is shown in Figure [l]^: ft rapidly decreases from 4Ao/e to about 
2Ao/e below He' = 40 mT, indicative of one lead going normal, whereas superconductivity 
in the second lead subsists up to He ^ 0.85 T, the junction being therefore of N-QD-S type 
at intermediate fields. 

For increasing coupling to the electrodes (i.e. higher conductance, assuming constant 
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FIG. 1: (a) SEM image of an aluminum nanogap obtained by electromigration of a constriction 
fabricated by angle evaporation. The scale bar is 300 nm. The nanometric gap created by the 
electromigration process is visible in the center, (b) Differential conductance map of device A as a 
function of gate and bias voltage, measured at zero magnetic field and T = 35 mK. (c) Extent of the 
non conducting bias voltage region at degeneracy gate voltage in sample A as a function of magnetic 
field, (d) Schematics of the level alignment for minimal bias voltage conditions (Vds — i2Ao/e) 
allowing for finite current flow. 

asymmetry and number of conducting channels) co-tunneling events can occur at low bias 
voltage and enhance the conductance via the Kondo effect at low temperature. The increase 
of the QD DOS by virtue of the Kondo effect originates from quasiparticles in the contacts 
lying close to the Fermi level, hopping on and off the dot. The Kondo effect and the 
suppression of the leads' density of states (DOS) at low energies caused by superconductivity, 
are therefore rather antagonistic processes. We illustrate the competition between Kondo 
correlations and superconductivity by two devices, B and C, with Kondo temperatures much 
larger and lower than Tc respectively. In the normal state of device B, a coulombian addition 
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energy Eadd > 80meV is observed to coexist with rather strong lead couphng, translating 
into a large Kondo resonance with = 14K (Figure [2]3, details in the Supplementary 
Information). In the superconducting state of the leads, superconductivity does not suppress 
the Kondo resonance and the superconducting features such as the coherence peaks at 
Yds = ±2Ao/e superimpose on top of it (Figure [2^). Such a behavior was predicted [T4l |24] 
and previously observed [25l [26j in other types of QD devices verifying Tk ^ Tc. The 
opposite limit Tk < Tc is illustrated by the more weakly coupled device C (Figure [S]^), 
which similarly displayed a Coulomb blockade and Kondo resonance pattern, but with a 
much lower Kondo temperature = 0.70 K. In the superconducting state, the Kondo 
peak at the Fermi level is destroyed [HI [25j. Interestingly, inside this CD, corresponding 
to an odd occupation number, i.e. to a local magnetic moment no longer screened by the 
lead electrons, the differential conductance is not maximum at Vds = ±2A/e. Further, a 
sign inversion of the differential conductance is observed, which is most prominent close 
to the degeneracy points. This negative differential conductance, as already reported in 
larger superconducting QD systems [23l [271 EH], results from Fermi's golden rule applied 
to the superconducting leads' DOS, which are convoluted with the delta function-like DOS 
of the QD. The gradual reduction of the superconducting gap under applied magnetic field 
reestablishes the Kondo resonance, which is further Zeeman split above about 400 mT 
(Figure |2|l). 

Moving on to high coupling strengths F, we show that Josephson supercurrent can be 
observed in SMTs. In equilibrium, dissipationless transport accross a S-QD-S structure is 
achieved by co-tunneling of Cooper pairs. The Josephson current amplitude Ic therefore 
scales with F^ and is experimentally only accessible for relatively strong coupling to the 
leads. Device D illustrates the molecular junction behavior in this regime. In the normal 
state (Figure |3^), the differential conductance map exhibits a Kondo resonance for positive 
gate voltages on a ridge delimited by only faintly visible CD edges. In this regime the 
electrostatic on-site repulsion is largely overcome by the strong coupling to the leads. The 
maximum differential conductance Gn ~ 0.85x (2e^//i) observed here approaches the unitary 
limit, while the width of the Kondo resonance is associated to a large Tk = 80 K, as 
already reported in normal SMTs in the past [121 EHJ- In voltage biased measurements of 
the superconducting state, the differential conductance around the Fermi level displays the 
usual coherence peaks at Vds = ±2A/e with A ^ 120 /xeV (not shown, similar to figure [2^). 



5 



For supercurrent measurements however, a current bias is better adapted. Figure [3)3 shows 
Vds{I) measurements for different gate voltages, showing the gate dependence of the critical 
current. In further detail. Figure [s]^ shows a dV/dl differential resistance map for small 
bias currents. The ridge visible here defines a switching current Ig such that d^V/dP |/^= 0. 
For |/| < Is the differential resistance is strongly reduced below Rj. We interpret the 
residual resistance below Is as related to phase diffusion of the Josephson state along the 
'tilted washboard' potential [cite]. No hysteresis is observed here which indicates that the 
Josephson junction is in, or close to, the overdamped limit. The strong reduction of the 
switching current with respect to the ideal critical current value Io{Rj) = A/eRj (equal to 6 
to 20 nA here, depending on the gate voltage) is generic to small Josephson junctions, and 
is understood as an effect of finite temperature and environmental coupling (e.g. to phonons 
[15]). In particular the critical current dependence on the normal state conductance Gn is 
expected 0, ^ to follow 

IJIo = (1 - V^l - GNh/2e^f/' (1) 

in an unshunted junction hosting a single spin-degenerate conductance channel. Experimen- 
tally, this 3/2 power-law dependence was observed to hold in metallic superconducting single 
electron transistors ^30j and in carbon nanotube junctions [7j. The fit of the Is{Gn) data 
shown in figure |3]l to (1) yields Iq = 1.01 nA. It is remarkable that whereas the justification 
of (1) strictly holds for the underdamped regime [30], a very good agreement to experiment 
is still found in our device D. 

In conclusion, we have shown a full characterization of Ceo SMTs connected to super- 
conducting leads, covering all coupling regimes. The data demonstrate the coexistence and 
competition of the effects of Coulomb repulsion, Kondo correlations and superconductivity 
in SMTs for various coupling strengths to the leads. In the most strongly coupled devices, 
the unitary conductance limit is approached in the normal state and for the first time a 
gate modulated superconducting current flowing through a single Ceo molecule is observed. 
The observed transport properties are in good agreement with the expected behaviour of 
a quantum dot connected to superconducting leads. In contrast with previous experiments 
involving carbon nanotube junctions [7j the energy spectrum is fixed by the molecular edifice 
itself and not by the nanoengineered portion of nanotube linking the contacting electrode. 
This experiment paves the way for further studies targetting at phase-sensitive measure- 
ments and the interplay of superconductivity with magnetically active molecules such as 

6 



endofuUerenes [5]. 
METHODS 

The starting material for the devices are thin metahic wires with a local constriction about 100 
nm wide and 20 nm thick, defined by e-beam lithography and double angle evaporation through a 
suspended mask, on top of a local AI/AI2O3 backgate (Figure[l^), as described in [13j. Ceo diluted 
in toluene is drop-casted on the thoroughly cleaned devices which are then immediately introduced 
into the sample mount thermally anchored to the mixing chamber of a dilution refrigerator. At 
liquid helium temperatures, the bias voltage across the wires is ramped up until the electromigration 
process sets on. The device resistance, which has an initial value of about 100 fi, increases, often 
stepwise, indicating the breaking up of individual ID atomic conduction channels. As soon as 
the device resistance exceeds 12.9 kfi, the bias voltage is cancelled within microseconds. After 
electromigration, the junctions are further cooled down to sub-0.1 K temperatures. About 50 % 
of all electromigrated junctions display a tunneling transport behavior with a low bias resistance 
Rk < Rj < I Gfi. About 10 % among these further show a gate dependent signal, indicative of QD 
behavior. High bandwidth control of the electromigration process is known to strongly increase 
the yield for obtaining narrow tunneling junctions. Two types of superconducting electromigration 
junctions have been investigated : (i) all aluminum electromigration junctions and (ii) proximity 
superconducting Al/Au junctions. In the later, the tunnel junction is opened through a short 
(L < 150nm) and thin gold wire connected to nearby larger superconducting aluminum source and 
drain pads (for details see Supp. Info.). The method used here for nanogap formation requires 
to get a low impedance environment in order to keep the bandwith during electromigration high 
enough during conductance monitoring. This forbids the implementation of chip-resistors therefore 
preventing a full control of the electromagnetic environment. 
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FIG. 2: (a) and (b) Differential conductance as a function of bias voltage at constant gate in two 
samples, at T below 40 mK and in both the normal (black line) and the superconducting state 
(red line), (a) Sample B has = 14 K, and superconductivity superimposes on top of the Kondo 
resonance. The inset shows the normal state data over a larger bias voltage window, (b) Sample 
C has Tk = 0.7 K, and the onset of superconductivity suppresses the Kondo resonance. The 
tunnel coupling of the BCS-density of states to a discrete level leads here to a negative differential 
conductance close to Yds — ±2Ao/e. Above 65 nS, the vertical scale is logarithmic, (c) Differential 
conductance plot of sample B in the normal state (T = 40 mK, B = 400 mT). The plots shown 
in (a) are taken along the dotted line (the color code can be read from the vertical scale in (a)), 
(d) Detailed magnetic field dependence of the low energy differential conductance in sample C. 
The coherence peaks merge into the Kondo peak at low fields, which is further Zeeman split above 
Be ^ 0.4T following Ez = ±gAS fiB {B - Be). As AS = 1, we find here g = 1.84 ± 0.2 along the 
best hnear fit (dotted lines). For better contrast, the differential conductance color code is different 
in the low and high field regions (the color code can be read from the vertical scale in (b)). 
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FIG. 3: (a) Differential conductance map of device D as a function of gate and bias voltage in 
the normal state (B = 20 mT). The dotted hues emphasize the CD edges, (b) Vds{I) in the 
superconducting state for = -3, -1.5, 0, 1.5, 3 V from black to green. No hysteresis is observed 
here, and the differential resistance does not go to zero below the switching current, but to some 
saturation value in the 5-10 kfi range, (c) Differential resistance map in the superconducting state 
(B = T) for small bias currents modulated by Ij^od = 18 pA at fo = 86 Hz. The green line is 
the zero bias dirrerential conductance in the normal state, (d) Switching current Is normalized to 
Iq = 1.01 nA as a function of normal state conductance (bullets), and fit to equation (2) (line). 
The only adjustable parameter is Iq. The data points deviating from the fit at high conductance 
correspond to measurements taken at ^ 4 V, where a significant tunnel leakage current from 
the gate electrode contaminated the measurements. 
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